Gap structure in noncentrosymmetric superconductors 
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Gap structure in noncentrosymmetric superconductors with spin-orbit band splitting is studied 
using a microscopic model of pairing mediated by phonons and/or spin fluctuations. The general 
form of pairing interaction in the band representation is derived, which includes both the intraband 
and interband pairing terms. In the case of isotropic interaction (in particular, for a BCS-contact 
interaction) , the interband pairing terms vanish identically at any magnitude of the band splitting. 
The effects of pairing interaction anisotropy are analyzed in detail for a metal of cubic symmetry with 
strong spin-orbit coupling. It is shown that if phonons are dominant then the gaps in two bands are 
isotropic, nodeless, and have in general different amplitudes. Applications to the ^(Pdi-^Pt^aB 
family of noncentrosymmetric superconductors are discussed. 
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I. INTRODUCTION 

Superconducting materials without inversion symme- 
try have recently become a subject of considerable in- 
terest, both experimental and theoretical. Starting 
from CePtaSi (Ref. 1), the list of noncentrosymmet- 
ric superconductors has grown to include Ulr (Ref. 2), 
CeRhSi 3 (Ref. S), CeIrSi 3 (Ref. !), Y 2 C 3 (Ref. 0), 
Li 2 (Pdi_ x ,Pt a; )3B (Ref. i), KOs 2 6 (Ref. 0), and other 
compounds. In most cases, the fundamental questions 
about the gap symmetry and the pairing mechanism re- 
main unresolved. 

The spin-orbit (SO) coupling of electrons with a non- 
centrosymmetric crystal lattice lifts spin degeneracy of 
the electron energy bands almost everywhere, which has 
important consequences for superconductivity: In the 
limit of strong SO coupling, the Cooper pairing between 
the electrons with opposite momenta occurs only if they 
are from the same nondegenerate band. This scenario 
is realized in CePtsSi, where the SO band splitting ex- 
ceeds the critical temperature by orders of magnitude. 8 
The same is likely to be the case in other materials, for 
instance in Li2(Pdi_ :E ,Pt, E )3B, see Ref. 0. 

The pairing interaction between electrons is most nat- 
urally introduced using the exact band states ^ 10 ' 11 ? 12 
which take into account all the effects of the crystal lat- 
tice potential and the SO coupling, see Sec. |TTJ In the 
strong SO coupling limit, the order parameter is repre- 
sented by a set of complex functions, one for each band, 
which makes the theory of noncentrosymmetric super- 
conductors similar to that of usual multiband supercon- 
ductors, see Ref. _13j. An alternative approach based on 
the representation of the pairing interaction in terms of 
the pure spinor states unaffected by the SO coupling was 
developed in Refs. HJ15| . 

In a phenomenological multiband pairing Hamiltonian, 
the relative strength of pairing in different bands can be 
arbitrary. In this article we go beyond the phenomeno- 
logical description and study the gap structure in noncen- 
trosymmetric superconductors under some fairly general 



assumptions about the microscopic mechanism of pair- 
ing. Specifically, we consider the interaction mediated by 
bosonic excitations (phonons and/or spin fluctuations). 
Starting with a microscopic expression for a momentum 
and frequency dependent pairing interaction, we derive 
the general form of the pairing interaction in the band 
representation, which contains both the intraband and 
interband pairing terms. The latter is shown to vanish 
identically in the case of isotropic BCS-contact interac- 
tion for any magnitude of the SO band splitting, see Sec. 
IIII1 In general, the interband pairing is absent only in 
the limit of large band splitting, see Sec. IIV1 

In Sec. El we present a detailed analysis of the possi- 
ble gap structures in noncentrosymmetric superconduc- 
tors of cubic symmetry, in a model which includes both 
the phonon and spin-fluctuation mediated interactions. 
The Conclusion contains a discussion of our results in 
the context of ^(Pdi-^Pt^B experiments. 



II. BASIC DEFINITIONS 

The Hamiltonian of noninteracting electrons in a non- 
centrosymmetric crystal has the following form: 

Hq = ^[e (fc)<5 Q/3 +7(fc)er Q/3 ]<4, Q a fc/ 3 

k 



(i) 



fe A=± 



where a, (3 —],[ are spin indices, a are the Pauli ma- 
trices, £\(fe) = eo(fe) + Aj'y(fe) j are the band dispersion 
functions, and the sum over fc is restricted to the first 
Brillouin zone. In Eq. ([T]) and everywhere below, sum- 
mation over repeated spin indices is implied, while sum- 
mation over the band indices is always shown explicitly. 
The SO coupling of electrons with the crystal lattice 
is described by the pseudovector ~f(k), which satisfies 
7(— fe) = — 'y(fe) and {g~f){g~ l k) = 7(fe), where g is any 
operation from the point group G of the crystal, see the 
examples below. 



The Hamiltonian in the first line of Eq. |T]) is diago- 
nalized by the following transformation: 



a-ka = J^ u a \(k)c k \, 

A=± 



(2) 



with the coefficients 




(3) 



V2|7l(l7l+A 7z ) 



forming a unitary matrix u(k). The Fermi surfaces de- 
fined by the equations £±(fe) = are split, except for 
the points or lines where 'y(fc) = 0. The band dispersion 
functions £\(k) are invariant with respect to all opera- 
tions from G, and also even in k due to time reversal sym- 
metry: the states |fe, A) and K\k, A) belong to k and — k, 
respectively, and have the same energy. Here K — i^Ko 
is the time reversal operation, and Kq is the complex con- 
jugation. One can write K\k, A) = t\(k)\ — k, A), where 
t\(k) = —t\(—k) is a nontrivial phase factor . 10 ! 11 For the 
eigenstates defined by expressions ([3]) we obtain: 



tx(k) = A 



lx{k) -ij y (k) 



ll(k) 



l 2 y{k) 



(4) 



The momentum dependence of the SO coupling is 
determined by the crystal symmetry. For the cubic 
group G = O, which describes the point symmetry of 
Li2(Pdi_ a; ,Pt a ;)3B, the simplest form compatible with 
the symmetry requirements is 



j(k) = 7 fc, 



(5) 



where 70 is a constant. For the point groups contain- 
ing improper elements, i.e. reflections and rotation- 
reflections, expressions become more complicated. In the 
case of the full tetrahedral group G = T^, which is rele- 
vant for Y2C3 and possibly KOS2O6, one has 



7 (fc) = l0 [k x {k 2 -kl)x+k y {kl-k 2 x )y+k z {k 2 x - 



e.) 



(6) 



This is also known as the Dresselhaus interaction,— and 
was originally proposed to describe the SO coupling in 
bulk semiconductors of zinc-blende structure. For the 
tetragonal group G = C4,,, which is relevant for CePtaSi, 
CeRhSi3 and CeIrSi3, the SO coupling is given by 

7(fe) = ~/±(kyX - k x y) + ^\\k x k y k z (kl - k 2 )z. (7) 

In the purely two-dimensional case, setting 7|| = one 
recovers the Rashba interaction^ which is often used to 
describe the effects of the absence of mirror symmetry in 
semiconductor quantum wells. 

Now let us take into account an attractive interaction 
between electrons in the Cooper channel, using the basis 
of the exact eigenstates of the noninteracting problem. 



The most general form of the interaction Hamiltonian 
the band representation is 



Hi, 



V Yl X! y A 1 A 2 A 3 A 4 (fc,fe';g) 



2V 



kk'q Ai, 2 ,3,. 



xc k+ q M c -k,X2 C - k '^ Ck '+i^- ( 8 ) 

We assume that the q-dependence of the pairing interac- 
tion is neglected (see the next section). The terms with 
Ai = A2 and A3 = A4 describe intraband pairing and 
the scattering of the Cooper pairs from one band to the 
other, while the remaining terms describe pairing of elec- 
trons from different bands. The above Hamiltonian can 
be considerably simplified in the absence of the intcrband 
pairing, which is the case if the SO splitting of the bands, 
Eso: is large compared with all energy scales associated 
with superconductivity. Since the pairing interaction is 
effective only inside the shells of width uj c (the cutoff en- 
ergy) in the vicinity of the Fermi surfaces, one can set 
Ai = A2 = A and A3 = A4 = A', and obtain: 



Hi-nt = ^EE^'( fe ' fc ') 



kk'q XX' 

XC fc+g,A ct -fc,A c -fc',A'Cfc'+q,A', (9) 



where 



V xx , (fc, k') = t x (k)t* x , (fc')VAA' (k, k') 



(10) 



The pairing amplitudes Vaa' are even in both k and 
k' (due to the anticommutation of fermionic operators) 
and also invariant under the point group operations: 
Vxx'(g- 1 k,g- 1 k') = Vxx'(k,k')^. 

In the case of large SO band splitting, the order pa- 
rameter has only intraband components. It is uniform 
(in the absence of external fields) and can be represented 
in the form A^(fc) = tx(k)Ax{k). The gap functions Aa 
transform according to one of the even irreducible rep- 
resentations of the point group and satisfy the following 
equations: 

AA(AO = -T£W|^AA<(fc,fc') 



n A' 



Av(fe') 



col + e y (k') + \Ay(k'W 



(11) 



The expression on the right-hand side converges due to 
the energy cutoff at uj c . 



III. BCS MODEL 

Let us calculate the pairing amplitudes and the gap 
functions in a simple BCS-like model, in which the at- 
tractive interaction is both instantaneous in time and 



local in space: 



Hint = -V / d 3 r^{(r)Vj(r)^(r)^ T W 



V 



- / d 3 r(ier 2 ) Q/ 3(icr 2 )^ 



X ^ ( r )^,T ^ ^ ^^ ^ ^ 



(12) 



where V > 0. Using the band representation of the field 
operators, 



v>, 



'a(r) = — =^2u aX (k)e lkr c kx , 

^ k,X 



(13) 



we obtain the pairing Hamiltonian in the form ([8]) with 
Vx 1 x 2 x 3 xA k , k ') = - — {ia 2 ) a p(i(r2)\s 

Here we neglected the difference between u a \(±k + 
q) and u a \(±k), which is 0(q/kp). In conventional 
centrosymmetric superconductors, we have q/kp ~ 
(ffcp) -1 -C 1 (£ is the correlation length). In the non- 
centrosymmetric case, the above estimate might not work 
and the q-dependence of the pairing interaction might be 
more important, leading, for instance, to the Lifshitz in- 
variants in the free energj&i^ and a spatial modulation 
of the order parameter even in the absence of external 
fields. We leave this issues to a separate publication. 
Using the identities 



u a x(-k) = t x (k)(ia 2 )ai3Up\{k), 



(14) 



and also the unitarity of the matrix u(k), we obtain for 
the pairing potential: 

Vx t x 2 x 3 x 4 (k,k') = --^tx 2 (k)tl 3 (k')Sx 1 x 2 S\ 3 \ i . (15) 

Therefore, interband pairing is absent in the BCS model 
for any strength of the SO coupling. Comparing this 
expression with Eq. (fT0|) . one can see that both the 
intraband pairing and the pair scattering between the 
bands are characterized by the same coupling constant: 
Vxx'(k,k') = — V/2. The pairing symmetry is isotropic, 
and it follows from Eqs. pip that the gap functions are 
the same in both bands: A + (fc) = A_(fc) = r\. This 
is not surprising, since the local interaction (TT2]) cannot 
lead to any fc-dependence of the gaps. 

The critical temperature is given by T c — 
{2e c /ir)uj c e~ 1 / NFV ', where C ~ 0.577 is Euler's constant, 
N F = (N + + N-)/2, and N x is the Fermi-level density 
of states in the Ath band. Although this has the usual 
BCS form, the superconductivity is non-BCS, because 
the order parameter resides in two nondegenerate bands, 
with T c and r\ independent of the band splitting and the 
difference between N+ and 7V_ . One can show that both 
the critical temperature and the gap magnitude are not 
affected by isotropic scalar impurities^ 



IV. INTERACTION MEDIATED BY BOSONIC 
EXCITATIONS 

Now we investigate a more general model, in which 
the pairing is assumed to be due to the exchange of 
some bosonic excitations. We consider two types of ex- 
citations: scalar (phonons), which couple to the elec- 
tron density p(r) = ipt(. r )'4'a{ r ), and pseudovector (spin 
fluctuations), which couple to the electron spin density 
s(r) = ^ a {r)cr a pipi3{r). Using the standard functional- 
integral representation of the partition function of the 
system, we obtain the following term in the fermionic 
action describing an effective two-particle interaction be- 
tween electrons: 

Sint = -^- / dxdx' p(x)D(x ~ x')p(x') 

+^j- dxdx's l {x)V lj {x-x')s 1 (x'), (16) 

where x = (r, r) is a shorthand notation for the coordi- 
nates in real space and the Matsubara time, J dx(...) — 

J dr J Q dr(...), g p h and g s f are the coupling constants 
of electrons with phonons and spin fluctuations, while 
D(x — x') and T>ij(x — x') are the phonon and spin- 
fluctuation propagators respectively. The spin fluctua- 
tions can be associated either with the localized spins, if 
such are present in the system, or with the collective spin 
excitations of the itinerant electrons (paramagnons)/21 In 
the latter case, T>ij{x — x') can be expressed in terms of 
the electron dynamical spin susceptibility XijiQ^)- I n 
general, the interaction (|16[) is nonlocal both in space and 
time. The BCS-contact Hamiltonian (Tj"2"| is recovered 
when the spin fluctuations are neglected and g^ h D(r,r) 
is replaced by — VS(r)8(r). 

In the momentum-frequency representation, Eq. (TT| 
yields the following pairing action: 

Sim = ^ Yl [9 2 P h D ( k - k')5 a s8fa 



kk' q 



+9 2 sf T) i3{k~- k )v a5 v ;. 



xa a (k + q)ap(—k)a- ( (~k')a$(k' 



Q), 



(17) 



where tt = 0V is the space-time volume, d a (k) and a a (k) 
are Grassmann fields, k = (fe,w n ), q = {q, v m ), and 
Lu n = (2n + l)nT and v m — 2rrntT are the fermionic and 
bosonic Matsubara frequencies, respectively. We assume 
that the conditions of the Migdal theorem are fulfilled, 
and also neglect the frequency renormalization, which 
corresponds to the weak-coupling limit of the Eliashberg 
theory. The theory developed below should work, at least 
qualitatively, even for such materials as CePtaSi, in which 
strong electron correlations are responsible for a heavy- 
fermion behavior and the above assumptions might be 
inapplicable. 

The phonon propagator is real and even in both fre- 
quency and momentum: D(k — k') = D(k' — k), and can 



therefore be written as follows: 

D(k - k') = D 9 (k, k') + D u (k, k'), (18) 

where the first term on the right-hand side, D 9 (k 1 k') — 
[D(k - k') + D(k + fc')]/2, is even in both k and k', while 
the second term, D u (k, k 1 ) = [D{k - k') - D(k + k')]/2, 
is odd in both k and k' . 

The spin-fluctuation propagator satisfies T>ij (k — k') — 
T>ji(k' — k) and can be broken up into the symmetric 
and antisymmetric in ij parts. Representing the latter 
in terms of a dual vector R, we obtain 

V l0 {k - k') = D» (fc, k') + D£.(fc, k') 

+i eij iRi(k - k'), (19) 

where the first (second) term on the right-hand side 
is an even (odd) function of k and k' , while Ri{k — 
k') = —Ri(k' — k). The antisymmetric component of 
the spin-fluctuation propagator is associated with the 
Dzyaloshinskii-Moriya interaction^ It is absent in the 
centrosymmetric case, due to the additional symmetry 
V tJ (k-k') =V lj (k' -k). 

After some straightforward algebra (see Appendix), 
the action (jTTJ) takes the following form: 



Sint = ^j^Vapjsik^k') 



kk' q 

xa a (k + q)a,0(-k)a^(-k')as(k' + q), (20) 

where the pairing interaction is represented as a sum of 
the /c-even, fc-odd, and mixed-parity terms: V = V 9 + 
yu _|_ ym _ rp ne even contribution is 

V^ s (k,k') = v g (k,k')(ia 2 ) al3 (ia2)l s , (21) 

where 

v g {k,k') = l -[gl h D 9 {k,k') - g 2 sf trV 9 {k,k>)}. (22) 
The odd contribution is 

Vap-ysik,^) = v u ^ J (k 1 k')(i(j l (j2) a p{iOj(J2)\ s , (23) 



where 



1 



v Utij (k, k') = - [g 2 ph D u (k, k') + g 2 sf trD u (k, k')]S tJ 

-g 2 s} Vl{k,k'). (24) 
Finally, the mixed-parity contribution is 

v ap-ys( k , k ') = v m ,i(k,k'){iaia 2 ) a fi{i^2)\ s 

+V m ,i(k', k){ia2) a ii{i^i(T2)\s, (25) 

where 

v m ,i(k, k') = -^-[Ri{k - k') + R t (k + k')]. (26) 



The first term on the right-hand side of Eq. (|25|) is odd 
in k and even in fc', while the second term is even in k 
and odd in k' . 

We would like to note that expressions (|2~lj) , (|23p and 
(|25| have completely general form in the sense that they 
do not rely on our assumptions about boson-mediated 
interactions and exhaust all possible spin structures of 
the pairing amplitude. Under the point group oper- 
ations g, the coefficients v g , v u jj, and v m transform 
like a scalar, a second-rank tensor, and a pseudovec- 
tor, respectively, and satisfy the invariance conditions 
v g {g^ 1 k 1 uj n ;g~ 1 k' ',w n >) = v g (k, u) n ; k', u n '), etc. By 
analogy with the theory of superconductivity in cen- 
trosymmetric compounds, see e.g. Ref. 24, Eqs. (|2T|) . 
()23|) correspond to spin-singlet and spin-triplet pairing 
channels respectively, while Eq. (|23|) describes singlet- 
triplet mixing. The possibility of singlet-triplet mixing 
due to the Dzyaloshinskii-Moriya interaction in the static 
case was pointed out in Ref. 1231 

Next, we use Eqs. @ to transform the pairing ac- 
tion into the band representation. Using identities (|14[) . 
we obtain the transformation rules for the pair creation 
operators in the spin-singlet and spin-triplet channels: 

(ia2)af3a a (k + q)ap{-k) 



Al,2 



tA 2 (fc)(5A 1 A 2 CA 1 (fc + g)cA 2 (-fc), 



(z<TCT 2 )a/3a Q (fc + q)ap(-k) 

] t\ 2 (fc)T Al A 2 (fe)cAi (k + q)c\ 2 (-k), 



Al,2 



where 

fi(fc) = u){k)&iu{k). 
Inserting these in Eq. ([2H|) . we obtain: 



(27) 



Sint = 2J| E E 1/ A 1 A 2 A 3 A 4 (fc^') 
kk'q Ai :2i3 .4 

x c Al (fc + q)c x . 2 {-k)cx 3 (-k')c X4 (k' + q), (28) 



where 

n i A 2 A 3 A 4 (fc,fc') -iA 2 (fc)n 3 (fc , )VA 1 A 2 A3A 4 (fc,fc'), (29) 

and 

VA 1 A 2 A 3 A 4 (fc,&') = v g (k,k')Sx 1 x 2 8x 3 x 4 
+v u ^j(k,k l )T^\ 1 \ 2 (k)Tj.\ 3 \ i {k') 

+V m .i(k, k')r l: x 1 X 2 (k)S\ 3 \ i 

+v m ^(k' ,k)5x 1 \ 2 T i ^ 3 \ i { k ')- (3°) 

The pairing amplitudes satisfy the following symmetry 
properties: 

^A 2 A!A 3 A 4 (-fc,fc') = AiA 2 V\ 1 A 2 A 3 A 4 (fc,^')j 
UA 1 A 2 A 4 A 3 (fc,-fc') = A 3 A 4 VA 1 A 2 A 3 A 4 (fc,&')- 



To obtain these, we used the anticommutation of the 
Grassmaim fields in Eq. Q28p and also the expressions 
dH) for the phase factors in Eq. ([29]) . 

It follows from Eq. (f50|) that, in general, all possi- 
ble channels are present in the pairing interaction, in- 
cluding interband pairing. The latter is absent, for any 
magnitude of the SO band splitting, if the odd harmon- 
ics of the bosonic propagators are negligible, so that 
v u ,ij(k,k') = and v m> i(k,k') = 0. This happens, in 
particular, for a fully isotropic interaction, in which case 
v g (k,k') = v g (u) n ,u> n >). 

We are particularly interested in the limit of large 
SO band splitting, which is relevant for the majority 
of noncentrosymmetric superconducting materials. In 
this limit, we set Ai = A2 = A and A3 = A4 = A' in 
Eq. ([50)) (the case of arbitrary band splitting, with both 
intra- and interband components of the order parame- 
ter present, will be considered in a separate publication). 
Since t\\ = A7(fc), the pairing action becomes 



s mt = ± i Y,Y,t x (k)t* x ,(k')vMk,k>) 

kk' q W 

xc x (k + q)c x (-k)cx>(-k')c x >(k' + q), (31) 



where 

Vxx, (k, k') =v g {k,k') + \X'v u< ij (fc, WYfi {k)% {k') 

+Xv m (k, k')j(k) + X'v m {k', m(k'). (32) 

This expression, together with Eqs. (J22J, tfM|) and ([2"6]) . 
relates the amplitudes of the intraband pairing and the 
interband pair scattering to the bosonic excitation spec- 
tra. Note that V\y (k, k') is even in both k and kl . Treat- 
ing the interaction (|31|) in the mean-field approximation, 
see e.g. Ref. [2J, one introduces the order parameters 
A x (k) = t\(k)A\(k,w n ), where, due to the symmetry of 
the pairing amplitudes, A\(—k,—u) n ) = A\(k,w n ). 



A. Weak coupling model 

In order to make progress, we approximate the fre- 
quency dependence of the pairing amplitudes by an 
anisotropic "square-well" models 

Vxx'(k,k') = Vxx>{k,k')e(u c -\un\)e{uJc-\v n >\) (33) 

where 9(x) is the step function, uj c is the frequency 
cutoff, and V\\/{k,k') depend on the directions of k 
and k' near the corresponding Fermi surfaces. The ap- 
proximation (|33|) has been used both for conventional 
phononic pairing interaction, see Ref. |25| . and also for 
spin-fluctuation mediated interaction, see Ref. |26J. The 
"square-well" decomposition also holds for the gap func- 
tions: A\(k,uj n ) = A x (k)9(ui c — \u> n \), so that the energy 
of quasiparticle excitations in the Ath band is given by 



Ex(k) 



The gap functions satisfy Eqs. (jlip . in which the Mat- 
subara sum is cut off at lo c . 

The pairing amplitude given by the matrix Eq. (|32|) 
is invariant under all operations from the crystal point 
group G, therefore Vx\>(g~ 1 k 7 g~ 1 k') = Vxx'(k,k'). 
Therefore, the momentum dependence of each matrix el- 
ement can be represented as a sum of the products of 
the basis functions of irreducible representations of G. 
In general, the basis functions are different for each ma- 
trix element. Neglecting this complication the pairing 
amplitude can be factorized as follows: 



Vxx>(k,k') 



E^v 



E 

4=1 



<Mfc)^*-(fc') 



(35) 



where a labels the irreducible representations (of dimen- 
sionality d a ) of G, which correspond to pairing channels 
of different symmetry, with (f> a ,i{k) being the even basis 
functions»2i The coupling constants V xx , form a Hermi- 
tian matrix, which becomes real symmetric if the basis 
functions are real. Keeping only the irreducible repre- 
sentation r which corresponds to the maximum critical 
temperature, the gap functions take the form 



i-h 



K x {k) =Ew^(fc), 



(36) 



and 57a, i are the superconducting order parameter com- 
ponents in the Ath band. The basis functions are as- 
sumed to satisfy the following orthogonality conditions: 
(cj)*(k)(j)j(k)) x = 5ij, where the angular brackets denote 
the averaging over the Ath Fermi surface. 

Linearizing the gap equations (|11[) we obtain the fol- 
lowing expression for the critical temperature: 



T r = 



2r L 



-uj c e 



-1/9 



where 




e x (k) 



|Aa(*0I : 



(37) 



g + -g-+ (38) 



is the the effective coupling constant, and 

gxx> = V xx ,N x ,. (39) 

While the critical temperature is the same for all dr com- 
ponents of r] X , the gap structure in the superconducting 
state below T c , see Eq. (|3"6"1) . is determined by the non- 
linear terms in the free energy, which essentially depend 
on the symmetry of the dominant pairing channel. 

V. PAIRING SYMMETRY IN A CUBIC 
CRYSTAL 

In the case of isotropic pairing interaction, one can 



(34) write v g (k,k') = v g (u> n ,u>n>) = -V g 9(uj c - \u n \)Q{uj c 



|k>n'|) in the square- well approximation. In this way, one 
recovers the BCS model of Sec. IIII1 with V = 2V g and 
the same isotropic gaps in both bands. 

To illustrate the effects of the interaction anisotropy 
on the gap structure, let us consider the following exam- 
ple. In a cubic crystal with G = O, the SO coupling can 
be described by ~f(k) = 70A;. This model is applicable 
to the Li2(Pdi_a;,Pt x )3B family of noncentrosymmetric 
compounds. The attractive interaction in these materi- 
als is likely mediated by phonons r^i therefore we ne- 
glect spin fluctuations by setting g s j = in expressions 
flU}, (El and (EH). Then, v g {k,k') = {gf }h /2)D9(k,k'), 
v u ,ij(k,k') = (g 2 ph /2)D u (k,k')5 lJ , and v m> i{k,k') = 0. 
Using the square-well approximation, one has 

v g (k, k') = v g (k, k')6(u c - \ui n \)0(u c - \uj n ' I), 
v u ,ij(k,k') = v Uii3 (k,k')9(uj c - \u n \)0(w c ~ |w„/|), 

with the momentum dependence inherited from the 
phonon propagator. Assuming a spherical Fermi surface 
and keeping only the s- and p-harmonics in the phonon 
propagator, we obtain: 



Vg(k,k')=-V g , 

Vu,ij{k, K ) = —V U (KK )Oij, 

v m ,i(k,k') = 0, 



(40) 



where V g and V u are constants. Note that this interaction 
is the same as the one considered phenomenologically by 
Edelstein in Ref. [lj. From Eq. ([32]) we obtain the pairing 
amplitudes in the band representation as follows: 



V X x'(k,k') = -V g - XX'V u (kk') 



(41) 



The components of the symmetric tensor kikj transform 
according to the representation A1+E+F2, where A±, E, 
and Fi are respectively one-, two-, and three-dimensional 
irreducible representations of the cubic group O (the no- 
tations are the same as in Ref. (28|). Therefore there are 
three pairing channels in the expansion (|35[) . with the 
following basis functions and coupling constants: 



V& = v 9+ o^A'K, 



^(fc) = l; 



(42) 



V& = ^AA'K, 



xx 15 



4> E {k) oc (k 2 x +ujk 2 +u*k 



2 , , j>; 



k'+^k'+uki); 
<Pp 2 {k) oc ykyk z , k z k Xi k x ky), 



where w = exp(27ri/3). 

Since phonons typically lead to a local attraction and 
cannot give rise to a substantial fe-dependence of the in- 
teraction, we expect that the Ai pairing channel dom- 
inates. Then the gap functions in the two bands [Eqs. 
([36]) ] are isotropic: A^(fc) = rjx, and satisfy the equa- 
tions 

Vx = Y^ 9 xx' ttT Y^ ■ 

A' 



x/^l + vi 



(43) 



where gxx' = Vxy Nx> ■ The critical temperature is given 
by Eq. ([3T]) . The gap magnitudes are not necessarily 
equal: For instance, in the vicinity of T c we find the 
following expression for the gap variation between the 
bands: 



(44) 



v+ - V- 


g++- 


-9— 


_ - 2g^+ + VV 


V+ + V- 


9++ - 


-9— 


_ + 2g_+ + VV 



where T> = \/{g++ — g ) 2 + 4^ g |_. Assuming that 

N + — A_ is small and that V g ^> V u , we have 



V u N + - A_ 

6Vg N F 



(45) 



Thus the gaps are different only if an appreciable p-wave 
harmonic is present in the phonon-mediated interaction 
and the SO coupling is sufficiently strong to create a 
considerable difference between the densities of states in 
the two bands. 

The coupling strengths being the same in both bands 
is not a generic situation. In the spirit of the standard 
model of two-band superconductivity— it is possible that 
the coupling constants corresponding to the intraband 
pairing channels and the interband pair scattering are all 
different. To obtain this we consider a generalization of 
the model (|4"0"1) which includes, along with phonons, also a 
contribution from spin fluctuations. In the absence of de- 
tailed information about the phonon and spin-fluctuation 
spectra in real noncentrosymmetric materials, in partic- 
ular in Li2(Pdi_ a; ,Pt a ;)3B, we use the model which in- 
cludes only the lowest angular harmonics consistent with 
the symmetry requirements: 



v g (k,k') 



-V 9 : 



j (k,k') = -V u (kk')8 i 

i\K, K ) = V m Ki. 



V u kikj, 



(46) 



Here the coefficients V g and V u are, in general, different 
from those in the model ([40]) . In the band representation, 
the pairing amplitudes become 



Vxx'(k, k') = -V g - XX' \V u (kk') 2 + K' 

-(X + X')V n , 



(47) 



There are three pairing channels, corresponding to the 
A\, E and F2 representations, see Eqs. (J42I) . The cou- 
pling constants in the A\ channel now have the following 
form: 



1 



v xx' = Vg + o XX 'Vu + XXX + (A + X')V m 



(48) 



The gap functions are isotropic: A.\(fe) = ijx, where the 
T]xs are found from equations (|43]) . The difference from 
the previous case is that now i] + ^= r\_ even if the density 
of states variation between the bands is negligible, i.e. 
N + = N- = Np. Assuming that V m is smaller than 



the other constants (i.e. the singlet-triplet mixing due to 
the Dzyaloshinskii-Moriya interaction is weak) , we obtain 
from Eq. (JSJ that 



V m 



V g -V u /3-V> 



(49) 



near the critical temperature. 

Finally let us consider the case of p-wave interaction 
dominating, which leads to an anisotropic pairing of the 
Fi symmetry. This happens if V u is large enough, and the 
degeneracy between the F 2 and E channels is lifted, e.g., 
by the Fermi surface anisotropy. The order parameter 
has the following form: 



A\(k) = \(riik y k z + rj 2 k z k x + r] 3 k x k y 



(50) 



The symmetry of the gap, in particular the location of the 
nodes, depends on the relation between the components 
of rj. There are four stable states of a three-dimensional 
order parameter in a cubic crystal: 29 (i) r\ — 770(1, 0,0), 
with two lines of nodes at k z = and k y = 0; (ii) r\ — 
770(1, i,0), with a line of nodes at k z = 0, and also point 
nodes at k x — k y = 0; (hi) 77 = 770(1, 1,1), with two 
lines of nodes at the intersection of the planes k x + k y + 
k z = ±1 with the Fermi surface, and also point nodes 
at k x = k y = 0, k y = k z = 0, and k z — k x = 0; and 
(iv) 77 = 770(1, w, w 2 ), with point nodes at k x = k y = 0, 



Ky — Kg ~~~ U f K% — K>£ - — - U . drUCl i^x — ^11 — ^z 



For 



the hrst three states one would have Cv(T) ex T 2 at low 
temperatures^ while for the last one cy{T) ex T 3 . 

It is instructive to interpret our results using the spin 
representation of the order parameter: 



A Q| a(fc) = ip(k)(ia 2 )ai3 + d(k){iaa 2 ) a f: 



where 



1&(fc) = - 



A+(fc) + A_(fc) 



is the spin-singlet component, and 



d(k) 



A+(fe)- A_(fc) 



7(fe) 



(51) 



(52) 



(53) 



is the spin-triplet component) 12 ! 31 The relative strength 
of the triplet and singlet order parameters is controlled 
by the difference between 77+ and 77-: |d|/|7/;| — r , see Eq. 
(I44p . In agreement with Ref. |15L only the component of 
d(k) which is parallel to ~y(k) survives (is "protected") in 
the limit of large SO band splitting. However, in the case 
of a weakly anisotropic phonon-dominated interaction, it 
follows from expression (|4"5")) that the triplet component 
is negligibly small. In the opposite case, when the in- 
teraction is strongest in the p-wave channel, one obtains 
from Eq. (j42|) that ip(k) = 0, i.e. the pairing is purely 
triplet. 



VI. CONCLUSIONS AND DISCUSSION 



We have studied the pairing symmetry in noncen- 
trosymmetric superconductors with SO splitting of the 
electron bands. The pairing interaction is derived using 
a microscopic model which includes both phonons and 
spin fluctuations. The interband pairing is shown to be 
absent for any strength of the SO coupling, if the interac- 
tion anisotropy is negligible. We have analyzed possible 
gap structures in the strong SO coupling limit with only 
intraband pairing and interband pair scattering present, 
using a cubic system as an example. If phonons are dom- 
inant, then the superconducting gaps in both bands are 
isotropic and nodeless (barring accidental zeros of the 
basis function of the unity representation), but do not 
necessarily have the same magnitude. 

Let us discuss the application of our results to the non- 
centrosymmctric compounds Li2(Pdi_ a; ,Pt. I ;)3B, where X 
ranges from to 1 (Ref. a). The critical temperature 
varies from 7-8 K for x = to 2.2-2.8 K for x = 1. The 
electronic band structure also exhibits considerable vari- 
ation: The SO band splitting in Li2Pd3B is as large as 
30 meV, while in Li 2 Pt 3 B it reaches 200 meV (Ref. [§), 
which in both cases is much larger than T c . Due to the 
absence of strong correlation effects and magnetic order, 
these materials provide a convenient testing ground for 
theories of noncentrosymmetric superconductivity. Su- 
perconducting pairing in Li2PdaB is due to the exchange 
of phonons, and the monotonic, almost linear, depen- 
dence of T c on the doping level x (Ref. M) suggests that 
it remains phononic for all x from to 1.— "21 

Experimental data on the magnetic penetration 
depth,— the electronic specific heat,—, and the NMR 
characteristics,— all seem to agree that Li 2 Pd3B is a con- 
ventional BCS-like superconductor with no gap nodes. 
In contrast, the gap structure in Li2PtsB is still a sub- 
ject of intensive debates. While earlier experiments, see 
Refs. 32.33.34, suggested the presence of lines of nodes 
in the gap, the recent /xSR and specific heat data 3 ^ have 
found no evidence of those. Moreover, according to Ref. 
l35l . the whole Li2(Pdi_ x ,Pt :r )3B family of compounds are 
single-gap isotropic superconductors. This conclusion is 
consistent with our results, see Sec. [V] Indeed, assum- 
ing that the pairing interaction in Li2(Pdi_ x ,Pt a; )3B is 
phononic and therefore only weakly anisotropic for all x, 
we obtain that the A\ channel always dominates, giving 
rise to nodeless isotropic gaps of essentially equal mag- 
nitudes in both bands. In order to create a noticeable 
difference between the gap magnitudes, see Eq. (|45|) . 
the interaction anisotropy would have to be very strong: 
Since (N + — N-)/Np ~ Eso/ e F and varies from 0.03 in 
Li2Pd 3 B to 0.2 in Li2Pt 3 B, the strength of the p-wave 
harmonic must be at least an order of magnitude larger 
than that of the s-wave harmonic, which is highly un- 
likely for a phonon-mediated interaction. 
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APPENDIX A: DERIVATION OF EQS. (|2lll26|) 



Let us start from Eq. (|17p . in which we substitute 
expressions (fTg|) and (TIT))) : 



S mt = ^J2{Spk[D 9 (k,k') + D u (k,k')}5, 



aSOfij 



kk' q 



+g 2 sf [V! j (k,k')+'D? j (k,k') 



+ie lJ iRi(k-k')}a l aS aj 3j } 
xa a (k + q)a f 3(—k)a~ / (~k')as(k' + q). (Al) 

The g-dependence of the fermionic fields plays no role 
in the algebraic transformations below, hence we use a 
shorter expression on the right-hand side: 



S int -^7^( I9 + IU + n . 



(A2) 



where 



19 = \Y,[9lhD 9 {k,k')5 aS 5p 1 +gl f V%{k,k l )^ aS al 1 } 



kk' 



x [a a {k)ap{-k) - ap{k)a a {-k)] 
x [a 1 (—k')as{k') — as(—k')a 7 (k')], 

/U = \Y.[tihD u {k,k')5 aS 5 (J1 +glp^{kX)<s°U 



kk' 



x [a a (k)ap(—k) + ap(k)a a (-k)] 
x [a 7 (— k')as{k') + as(— k')a~,(k')], 

I" 1 = h^iWlf £[««(* - fc ') + -**(* + k')]a l aS a- 



3 

01 



kk' 



x[a a (k)ap(-k) + ap(k)a a (-k)} 
x [a 7 (— k')as(k') — as(—k')a 7 (k')] 

+ I^vi9 2 s f E^( fc - fc ') - R l( k + k ')\<s<r, 



3 

01 



kk' 



x [a a (k)ap(-k) - ap(k)a a (~k)} 
x [a^(-k')as(k') + a s {-k')a 1 {k 1 )]. 



The even in k combinations of the fermionic fields can be 
represented as follows: 



a a (k)ap(-k) - ap(k)a a (-k) 

= -{i^2) a p{icr2)^ v a^{k)a v {-k), 

a 7 (—k')as(k') — as(—k')a 1 (k') 

= -(ia 2 )- ( s(io-2)l (7 a P (-k')a< 7 (k'), (A3) 



while the odd combinations have the form 



a a (k)ap(~k) + ap(k)a a (—k) 

= (*cr i cr 2 )^(«cr i cr2)piyap(fc)a l/ (-fc), 
a~f(—k')as(k') + as(—k')a~ / (k') 

= (io-i<T 2 )~ f 5(icri<J2)i lT ap(-k')a <T (k'). (A4) 



Using the matrix identities 



d a 8Sp 1 (ia2y af3 (ia 2 )- l s = 2, 

(&i)as(o'])f3 1 (icr2) a 0(io'2)-fS — — %5ij, 

S a 5Sf3 1 (ia i a2) ,i 3 (io-ja2) 1 s = 2Sij, 

{&i)a&{&j)l3' 1 {iOmV2) cl ( i {iOn&2) 1 5 
{Ci)a8(cj)(3~/(i(TmCr2) al3 (i<T2)~/5 = 2ieij m , 

(ci)as(o'j)0j(icr2) a p(i(T m cr2)"/6 = -2«e 



(A5) 



ijm , 



we arrive at Eqs. (J2HI26I 
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